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^ ABSTRACT 
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^ We present an interferometric kinematic study of morphologically complex proto- 

stellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The 
^ strong asymmetric nature of most envelopes in our sample leads us to question the 

^ common interpretation of velocity gradients as rotation, given the possibility of pro- 

jection effects in the observed velocities. Several "idealized" sources with well-ordered 
^ velocity fields and envelope structures are now analyzed in more detail. We compare 

the interferometric data to position- velocity diagrams of kinematic models for spheri- 
^ cH cal rotating collapse and filamentary rotating collapse. For this purpose, we developed 

Qh a filamentary parametrization of the rotating collapse model to explore the effects of 

O geometric projection on the observed velocity structures. We find that most envelopes 

^ in our sample have PV structures that can be reproduced by an infalling filamentary 

Ci envelope projected at different angles within the plane of the sky. The infalling filament 

produces velocity shifts across the envelope that can mimic rotation, especially when 
^ viewed at single-dish resolutions and the axisymmetric rotating collapse model does 

not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect 
^ rotation, then the inferred centrifugal radii are quite large in most cases, indicating 

^ significant fragmentation potential or more likely another component to the line-center 

^ — 5 velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation 

O alone when envelopes are non-axisymmetric and that projected infall velocities likely 

^ dominate the velocity field on scales larger than 1000 AU. 
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1. Introduction 



The angular momentum of protostellar cores and envelopes plays an important role in star 
formation. To the extent that the rotation of the protostellar cloud is non-negligible, a significant 



fraction of the star's final mass must be accreted through a circumstellar disk (Zhu et al. 2010 



Kratter et al. 2010). In addition, collapse with large amounts of angular momenta can result 



in gravitational fragmentation into companion stars (Burkert & Bodenheimer 1993; Boss 1995 



and references therein). Finally, differing initial angular momenta of clouds are likely to produce 
a distribution of binary separations and differing initial disk configurations, with implications for 
planet formation. Thus, the importance of angular momentum in star and planet formation process 
illustrates the need to characterize its magnitude in collapsing protostellar cores. 

Observational attempts have been made to characterize the rotation of dense cores and proto- 
stellar envelopes on scales of 0.05 — 0.5 pc, principally using the dense gas tracers NH3 and N2H+ 
( Goodman et al.||l993 ; Caselli et al.||2002 ). Velocity gradients are clearly detected in most objects 
and observations with increasing spatial resolution tend to find larger velocity gradients on the 
same spatial scales due to less smoothing of the velocity field. On smaller scales, [Volgenau et al 
( |2006) and |Chen et al. (2007) examined the kinematic structures of protostellar envelopes down 
to ^1000 AU scales, only finding clear evidence of rotation in a few cases. The detected velocity 
gradients are often interpreted under the assumption of solid-body rotation, even in protostellar 
systems; this enables simple calculations of angular momentum ( Goodman et al.||1993 ; Chen et al. 
2007; Chiang et al.||2010 ; Tanner fc ArcepOll ). Several studies have also found the possible signa- 
tures of rotation on sub- 1000 AU scales from high- velocity molecular line wings in interferometric 
data dBrinch et aLl|2007t |Lee et al.||2009D . 



In addition to rotation, the velocity fields of protostellar envelopes will necessarily include infall 



motions. Initial models of collapsing protostellar envelopes assumed spherical geometry ( ^Larson 
1969; Shu|1977 ); later models included angular momentum, which leads to flattening of the envelope 



on small scales near the forming disk (Cassen &; Moosman 1981; Terebey et al. 1984). Infall is 



typically probed via single-dish spectra, and the blue-shifted asymmetry of an optically thick line 



profile (e.g. HCO^ HCN, and CS) is the classic signature of collapse (Walker et al. 1986; Zhou 
1992; [Zhou et al.||1993 ). The inside-out collapse model ( Shu||1977 ) is typically adopted to describe 
the velocity field and coupled with radiative transfer modeling to derive infall rates and infall radii 



dZhou et al.]|1993t [Ward-Thompson fc Buckley||200H [Narayanan et al.||2002D . The derived infah 
radii are often on the order of ^5000 AU, indicating that material is likely falling-in from relatively 
large radii. Inside the infall radius, the envelope is expected to exhibit a differential rotation curve, 
owing to the conservation of angular momentum, the outer core is often assumed to have solid- 
body rotation f Terebey et al. 1984). However, Ward-Thompson & Buckley (2001) noted that one 



of their models required inward motions of the envelope out to its boundary, indicating that inside- 
out collapse may not be realistic and an "outside-in" collapse may be more appropriate in some 
cases(e.g. Foster fc Chevalierl[1993 ). 
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Contrary to the assumptions of early models, simulations of star formation in turbulent clouds 
indicate that infall to protostars is generally complex, non-axisymmetric, and often filamentary (e.g. 



Bate et al. 2003; Smith et al. 2011). Indeed, observations at higher spatial resolution show that 



most pre and protostellar cores are elongated (Benson &; Myers||1989; Myers et al.||1991 ; Bacmann 



et al.||2000t |Stutz et al.||2009D and more hkely to be prolate ( |Ryden|[T996l |Myers||2005l ) . In the first 
paper of this series (Tobin et al. 2010|), hereafter Paper I, we used Spitzer IRAC images of a set 



of Class protostars to characterize dense envelope structures in extinction at 8/xm; the results 
showed that many systems have very flattened, filamentary, and often strongly non-axisymmetric 
density structures. 

The recognition of the complexity and non-axisymmetry of many protostellar envelopes poses 
substantial problems for any attempt to infer angular momenta and the interpretation of line profiles 
under the assumption of spherical collapse. Kinematic observations necessarily refer to integrations 
along the line of sight of structures with uncertain geometry, and because only one of the three 
velocity components is detected, there is substantial ambiguity in interpreting the observations as 



rotation or infall (e.g. Dib et al. 2010[ Smith et al. 2012 submitted). More specifically, infall (or 
outflow) could masquerade as rotation, depending upon the geometry and viewing aspect of the 
envelopes. 



In the second paper of this series Tobin et al. (2011), hereafter Paper II, we presented obser- 



vations of the Class envelopes of Paper I in the dense molecular tracers N2H+ and NH3. Using 
moment maps, we showed that the kinematics of these envelopes, like their density structures, were 



often complex and difficult to interpret. The interferometric data presented in Tobin et al. (2011) 



had a higher frequency of envelopes with velocity gradients that could be interpreted as rotation. 



as compared to Volgenau et al. (2006) and Chen et al. (2007); however, we were not convinced that 



our data were clearly tracing a rotation signature, even in ideal cases. 

In this paper, we proceed to a more detailed analysis of the kinematics for a subset of the 
sample from Papers I & II where the velocity fields appear to be rather well-ordered. Mindful of 
the interpretation difficulties in complex geometry, we develop a simple kinematic model for the 
collapse of a filament to compare with the observations. While significant ambiguities remain, we 
find evidence for rotation in some systems on small (< 1000 AU) scales, but infer that infall motions 
are likely to dominate the radial velocity structure on scales > 1000 AU. Our findings emphasize 
the importance of following the kinematics to the smallest scales to determine the angular momenta 
of infalling envelopes, and the need for further, more realistic simulations of the large-scale velocity 
fields in star-forming regions for comparison with observations. 

This paper is organized as follows: Section 2 briefly describes the models of filament and ax- 
isymmetric collapse. Section 3 describes the comparison of the models to the observed channel maps 
and PV diagrams, and Section 4 discusses our results. The observational details were described in 
Paper II and are omitted here for brevity. We will present data for five sources, listed in Table 1, 
four of which observed the N2H+ molecule (one with the Plateau de Bure Interferometer and three 
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with CARMA) and one source was observed in NH3 with the Very Large Arra}Q 

2. Kinematic Models 

Most previous work on the velocity structure of protostellar envelopes assume spherical or 
axisymmetry in the interpretation of kinematic data. However, in the absence of symmetry, the 
physical interpretation of velocity gradients is now non-trivial. The difficulties of interpreting 
velocity fields of protostellar envelopes seen in projection are schematically illustrated in Figure 1. 
For simplicity we plot only the limiting radial velocities expected from the emitting envelope. In 
the case of spherical symmetry (or axisymmetry in the appropriate plane), pure uniform rotation 
(solid-body) directly yields a simple linear velocity gradient across the envelope; the velocity spread 
is due to the superposition of velocity components at intermediate scales and disappears toward the 
envelope edge. Pure radial infall yields no systematic gradient across the envelope, but an increase 
in the velocity dispersion near the central gravitating mass. In general we expect both types of 
motions to be present, resulting in the skewed velocity caustics shown in the upper right. Note that 
a Keplerian rotation curve is distinct from that of solid-body rotation and appears exactly like the 
radial infall curve. 

Proceeding to a filamentary envelope, the simplest of non-axisymmetric structures, immedi- 
ately introduces a fundamental ambiguity, as both pure solid-body rotation and pure infall can yield 
similar velocity gradients from one side of the envelope to the other (bottom two panels). Including 
both types of motion makes the situation even more uncertain. As the observations clearly indicate 
non-axisymmetric or filamentary envelopes are common (Paper I), this poses a serious problem for 
interpretation. 

In principle, the emission in position-velocity (PV) diagrams, such as indicated schematically 
here, can provide some indication of true filamentary or non-axisymmetric structure. For instance, 
axisymmetric envelopes (especially those infalling) will tend to have larger velocity dispersions at 
a given position because material along the line of sight will have differing velocity vectors than in 
a highly linear structure. In addition, material falling-in close to the protostar will tend to yield 
both red- and blue-shifted emission at each position, whereas a filament will exhibit more of an 
abrupt change from red to blue across the position of the protostar (bottom right panels of Figure 
!)• 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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2.1. Analytic Model of Filamentary and Axisymmetric Collapse 



Given the complications introduced to the kinematics by non-axisymmetric systems, we make 
an initial step toward interpreting the observations by developing a simple kinematic model for 



comparison with position- velocity diagrams. We have modified the rotating collapse model (Ulrich 



1976; Cassen fc Moosman||1981, hereafter CMU model) which has been extensively used in studies 



of protostellar envelopes. The infalling gas in the CMU model is assumed to fall in from an initial 
spherical cloud rotating with constant angular velocity. With this assumption, the streamlines of 
the infalling gas do not intersect, and thus the (supersonic) motion can be considered in the limit 
of ballistic trajectories around a central gravitating mass. (This neglects the self-gravity of the 



envelope; but Terebey et al. (1984) showed that the CMU model could be taken as the inner region 



of a more general model of self-gravitating spherical cloud collapse.) 

By neglecting the self-gravity of the infalling gas, we can adopt the axisymmetric velocity field 
of the CMU solution. We create a filament by modifying the density distribution such that all the 
infalling gas is confined within specified streamlines spanning A(/)=30° in azimuth on opposite sides 
of the envelope (mirror symmetry). This filamentary geometry enables us to project the density 
distribution of the envelope to different orientations in the sky in order to examine the effects of 
projection on the kinematic structure. The height of the envelope in the z-direction is taken to be 
comparable to the thickness of the filament as viewed from top-down, but does not play a crucial 
role in the observed kinematic structure. 

The streamlines comprising the filament correspond to limits in velocity space for emission at a 
given position, producing curves which can be compared with observed position-velocity diagrams, 
a similar concept to that of using "caustics" to analyze position-redshift diagrams of galaxy clusters 
(e.g. Regos &; Geller 1989). The velocity components of the infalhng, rotating gas are given by 
( |Uhich|[T976l ) 
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describing parabolic motion around a central point mass. The angle is the angle between the 
orbital plane and the rotation axis, while 9 is the angle from the rotation axis to the particle. For 



* Ulrich (1976) and Chevalier (1983) both have sin in the denominator, rather than sin^^ as written in Terebey 
et al. (1984) and Hartmann (2009). The correct term is sin^^; this can be verified by ensuring that the squares of 
the three velocity components add to 2GM/R. 
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simplicity, the velocities are only considered for motion nearly in the equatorial plane (^o ^ 89°), 
making vq ^ 0. Note that the term cos 9/ cosOq, appearing in both Vr and V(f) will go to zero as 

0. Furthermore, most envelopes in this study are rather filamentary and we will be comparing 
with PV diagrams taken from equatorial regions, making vq contributions negligible. A fiducial 
central object mass of 0.5 Mq is assumed in these models, intended to be a typical protostellar 
mass in the absence of any real constraints. However, the envelope masses within 10,000 AU are 
^1 Mq for the objects that will be considered in this work (Table 1); this is not negligible with 
respect to the assumed central mass and we will discuss its possible effects on the kinematics in 
Section 4. 

Other than mass, the only free parameter of the CMU velocity field is the initial angular 
velocity of the material falling in at the current time. Material falling in from an initial direction 
given by the angle from the rotation axis lands on the (flat) disk at a radius r, given by 

JL = (4) 

Rc 1 -cos6'/cos(9o ^ ' 

where Rc is the centrifugal radius, 

i?o is the radius at which the material fell-in from with angular velocity ft. The cos 9/ cosOq terms 
in Equations 1-3 can then be rewritten in terms of r/Rc, relating the CMU rotation velocity {V(p) 
to the initial angular velocity (fi) of the infalling shell. 

To calculate the velocities that would be observed in a PV plot across the equatorial plane, 
the velocities in the r, 0, and 6 directions must be projected along a line of sight. For simplicity, 
the line of sight is defined to be the y-axis at z = 0, which makes 

vios = i;r(sinzsin$o) + veicosisin^o) + 'u^(cos$o) (6) 

where $o is projection angle of the envelope within the line of sight, measured counterclockwise 
from the x-axis. $o is unrelated to (f) which is the azimuthal angle of infalling material and $o 
only affects the observed velocity distribution in a non-axisymmetric system. The angle i is the 
inclination of the rotation axis of the system with respect to the plane of the sky; this adds 
an additional projection term to the observed velocity distribution. However, we only consider 
velocities calculated at inclination of 90°, meaning that the rotation axis of the model would be in 
the plane of the sky. This simplification is made to limit free parameters and because most of the 
observed systems in Paper II have inclinations of z > 60°, making this projection effect minimal. 

We neglect radiative transfer and chemistry for this initial exploration; this is reasonable 
because we will compare the models to optically thin molecular line kinematics that trace cold 
(T^IOK) molecular gas at the densities expected for protostellar envelopes (n > 10^ cm~^, i.e. 
NH3 and N2H+). These molecules are present with detectable emission from large-scales (^ 0.1 
pc) down to radii as smah as -1000 AU dCaselh et aLl|2002t [Benson fc Myers||1989t |Lee et al.||2004 
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Paper II); however, the abundances of these molecules are not constant with radius. Our models 
use the simplest possible non-axisymmetric envelope structure to describe the high density material 
and optically thin emission line kinematics: an m = 2 single filament. The use of optically thick 
lines to probe the envelope kinematics would need to consider lower density surrounding material 
that could contain even higher-order structure resulting in more complex line profiles. Comparison 
to the numerical models of Smith et al. (2011 ) with radiative transfer will be the subject of a future 
paper. 



2.2. Predicted Kinematic Structure 

The left panel of Figure [2] shows the infall streamlines on large scales for the axisymmetric case, 
and the middle panels show a zoom-in on the inner 1000 AU. The infall streamlines are computed 
every 15 degrees in the (j) direction to sample the entire envelope. The projected velocities versus 
X-axis position are used to construct simple PV plots from the individual envelope streamlines, 
shown in the right panel of Figure [2j Note that we have assumed that the envelope is infalling 
from 10,000 AU, with no static outer core. However, we know that there is substantial mass on 
larger scales for all these objects measured via 8/xm extinction, continuing to increase from 0.05 pc 
to 0.15 pc in radius (Paper I). The masses compared with the observed linewidths indicate that 
the envelopes are gravitationally bound (Paper II); therefore, we expect that material from larger 
scales will be falling in, but under the influence of greater enclosed mass. Thus, the large-scale 
velocities from the CMU velocity field are lower limits. 

To approximate a filamentary envelope in the same manner as the axisymmetric envelope, we 
only calculate three streamlines on either side of the protostar, spanning A(/)=30° in azimuth at 
scales larger than Rq- The geometry of the infall streamlines is shown in the top left and middle 
panels of Figure [3] The PV plots are then constructed in the same manner as the axisymmetric 
case and shown in the right panels of Figure [3] We have assumed that matter is infalling from 
10,000 AU, the same as the axisymmetric case. Notice that the region of velocity space enclosed by 
the filament model is smaller than the axisymmetric case; the velocity streamlines plotted in the 
filament model represent a subset of those plotted for the axisymmetric model. Then the effects 
of viewing the filamentary envelope at $o = ±30° projection angles are shown the middle and 
bottom rows of Figure [3j $o specifically refers to the projection angle of the middle streamline on 
either side of the envelope in Figure [3j Variations of $o in non-axisymmetric systems can clearly 
introduce apparent velocity shifts and gradients that are not present in axisymmetric geometries. 

Substantially different velocity structures are evident in the filamentary PV plots, especially 
when viewed with different projection angles within the plane of the sky, as compared to axisym- 
metric collapse in Figure [2| Notably the velocity shift across the filamentary envelopes on scales 
> 1000 AU is not due to rotation velocity; rather it is due to the infall velocity being projected 
along the line of sight. The schematic plots of envelope kinematics in Figure 1 and the PV plots 
derived from the CMU model for axisymmetric and filamentary envelopes can he used a benchmarks 
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to compare with our observational data. This will potentially enable us to distinguish between ax- 
isymmetric and non-axisymmetric inf ailing envelopes , as well as velocity gradients that are due to 
rotation or projected infall. 

3. Observational and Model Results 

We have selected five protostellar systems from Paper II to compare with the kinematic models, 
all having relatively uncomplicated morphological structure and/or ordered kinematic structure. 
These features are important because the interpretation of the velocity structures is complicated 
by the often complex envelope structures (Papers I & II). We examine channel maps and position- 
velocity (PV) diagrams for each protostellar system to better understand the dynamical processes 
behind the observed kinematic structure. We compare the models and data using PV plots of 
the simplistic axisymmetric or filamentary collapse models described in the previous section. The 
final model parameters are listed in Table 2. These models do not provide unique fits, given the 
geometric uncertainties, but they illustrate the basic features and issues involved in determining the 
processes giving rise to envelope kinematics. Also note that we have not combined zero-spacing data 
with our interferometric data resulting in the large-scale emission being filtered out and possibly 
influencing our kinematic data. However, we show in Paper II that the velocity fields observed in 
the single-dish and interferometric data are consistent at large-scales, indicating that the spatial 
filtering by the interferometer is not biasing our kinematic results. 

3.1. L1157 

L1157 is located in Cepheus at a distance of ^300 pc and is an example of a rather well-ordered, 
highly flattened envelope. The flattened envelope is observed in Sfim extinction, as well as N2H+ 
and NH3 emission. The N2H+ integrated intensity map is shown in the upper left panel of Figure [ij 
The envelope mass was found to be at least ^0.86 Mq from 8 /xm extinction (Paper I) and possibly 
an additional ^0.7 Mq inferred from the N2H^ emission (Paper II). We have further shown in 
Paper I that the flattening reflects a filamentary structure in three-dimensions rather than a sheet 
extended within the plane of the sky. The scale heights of the short envelope axis were ^1600 AU 
across all radii (except the inner iC), consistent with a hydrostatic filament; a sheet would be a 
factor of three more extended in the vertical direction. Thus, L1157 may be well-described by an 
axisymmetric filament, one of the ideal cases shown in Figure [TJ 

The channel maps in Figure |4] show several kinematic features in this system. In the top and 
bottom rows, there is higher-velocity N2H^ emission north and south of the protostar extended 
along the outflow. This emission was attributed to outflow interactions in Paper II. In the middle 
panels, there are emission peaks from the inner envelope on ±5^' scales present in all the channels. 
However, on scales > 7'' we see that the east side of the envelope comes into view at blue-shifted 
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velocities and the west side at red-shifted velocities. This traces a clear large-scale velocity gradient 
that may be attributed to rotation or projected infall. Note that the extension of N2H+ emission 
southeast of the protostar traces emission from the outflow cavity wall and also appears to reflect 
outflow entrainment. 

The PV plot in Figure [5] also shows high- velocity emission due to the outflow-envelope inter- 
action extending toward red-shifted velocities on both sides of the envelope, coincident with the 
N2H^ emission peaks. Thus, any indication of high-velocity emission due to rotation or infall at 



small-scales is masked by the outflow effects. However; |Gueth et al. (1997) appears to have detected 
a smaller-scale velocity gradient in C^^O emission On scales >T' in Figure [s) the envelope velocities 
are nearly constant with a systematic ^0.1 km s~^ velocity shift between the east and west sides 
of the envelope, as seen in the channel maps. The larger-scale features, with their nearly constant 
velocity, are quite similar to the filament model shown in the bottom right panel of Figure [3| At 
the same time, the velocity structure is dissimilar to uniform rotation and the axisymmetric infall 
model (c.f. [Chiang et al. 2010|). A filament model has been overlaid on the data in the PV plot 



in Figure [5| showing that the fit is reasonable, except in the inner envelope where the kinematic 
structure is not accurately probed. The model is projected within the plane of the sky by $0 — 
15° and has a centrifugal radius of 100 AU (Table 2). 



3.2. L1165 

LI 165 is also located in the Cepheus region and the protostar is forming within a larger- 
scale filamentary structure (Paper I). The N2H+ emission on small scales, shown in the upper left 
panel of Figure [6} traces a structure extended normal to the outflow, peaking just southeast of the 
protostar. A lower limit on the envelope mass is measured to be ^1.1 Mq from 8/xm extinction in 
Paper I, with perhaps at least an additional ^0.4 Mq inferred from N2H^ emission in Paper II. 
The N2H^ channel maps for LI 165 are shown in Figure |6j indicating a clear velocity gradient across 
the protostar (normal to the outflow direction) . The velocity on the northwest side of the protostar 
is blue-shifted and nearly constant. The emission then becomes red-shifted on the southeast side 
of the envelope, with the highest velocity emission being red-shifted, near the protostar. There is 
evidence of some higher-velocity blue-shifted emission adjacent to the protostar on the northwest 
side in the channel maps, but it is not as definitive as the red-shifted emission. 

The PV diagram in Figure [t] shows that the N2H^ emission northwest of the protostar has a 
very narrow linewidth and a roughly constant velocity; the broadest linewidth is southeast of the 
protostar, near the N2H^ peak. Southeast of the peak linewidth, the line rapidly becomes narrow 
again, while becoming more blue-shifted. The channel maps show that the N2H+ kinematics are 
unlikely to be outflow related, due to their location orthogonal to the outflow. The PV structure in 
Figure [7| is most similar to the filamentary infall model; its velocities are overlaid in Figure [7| The 
model is able to approximate the observed velocity structure, though there are discrepancies; for 
LI 165 the model is projected by $o=15° within the plane of the sky and has a centrifugal radius 
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of 10 AU (Table 2). The small centrifugal radius of the model implies that the data are consistent 
with little rotation on scales probed by N2H+. 

As Figures [g] and [t] show, N2H+ is only able to trace emission to within R^1200 AU of the 
protostar. A special feature found in LI 165 is that HCO^ (J = 1 ^ 0) emission is able to trace 
small-scale kinematics with higher velocity in the inner envelope, shown in Figure [8j The red 
and blue-shifted components are extended normal to the outflow at radii of 3^' (600 AU) from the 
protostar. This emission lies inside the broadest N2H^ emission, suggesting that it traces smaller- 
scale kinematic structure, in agreement with the small-scale model velocities which are overlaid 
in Figure |8| Furthermore, location of the blue and red-shifted high- velocity emission on opposing 
sides of the envelope with opposite velocities, along with the small velocity width a large-scale in 
the N2H+ emission, is consistent with what is expected for a filamentary envelope (Figure [s]). A 
similar high- velocity HCO^ feature was found in RN043 (Paper II). 



3.3. CB230 

CB230 was classified as a "one-sided" envelope in Paper I, given its strongly asymmetric 
distribution of 8/xm extinction; this was reflected on large-scales in the single-dish N2H+ map 
shown in Paper II. We measure a lower limit on the envelope mass within 0.05 pc of ^1.1 Mq 
in Paper I from 8/xm extinction, but there is clearly mass on small-scales that we do not probe. 
The smaller-scale emission probed by NH3, shown in the upper left panel of Figure [9) is slightly 
asymmetric, but well-ordered as a whole. The mass inferred from the small-scale NH3 emission in 
Paper II was ^4.8 M©, but highly uncertain. There is also a depression in the NH3 emission in 
the inner envelope, coincident with the protostar. This effect is attributed to destruction of NH3 
in the inner envelope (Paper II). The velocity structure of the NH3 emission is shown by channel 
maps in Figure |9| Because the main NH3 (1,1) lines are a pair separated by 0.2 km s~^, the last 
three panels in the middle row are blends of blue and red-shifted emission, while the rest of the 
plots reflect mostly unblended emission. The channel maps show that there is a clear, well-ordered 
velocity gradient across the envelope from east to west, normal to the outflow. 

The PV diagram in Figure [TOj in contrast to the channel maps, is derived from NH3 satellite 
emission lines that are separated by ^0.4 km s~^. This makes the emission appear more distinct 
than if the PV diagram were plotted using the main NH3 lines. The PV plot and channel maps 
show that the emission on the east and west sides of the envelope has a rather constant velocity and 
there is an abrupt shift from blue to red-shifted emission starting at the location of the protostar. 

While a measure of the small-scale velocities is missing, the emission on scales >1000 AU is 
most consistent with the infalling filament, due to the constant velocity emission out to large-scales 
and the relatively narrow linewidth. The constant velocity emission is dissimilar from what would 
be expected for uniform rotation, which would have linearly increasing velocity out to large-scales. 
The filament model is overlaid on the data in Figure [TO) having projected by $0 = 20° and has Rc 



= 10 AU. The model does not exactly describe the data, but substantial mass in the inner envelope 
implied by the NH3 data could make the velocity shift slightly larger on the scales probed. 



3.4. IRAS 16253-2429 

IRAS 16253-2429 is the most nearby object in our sample, located in Ophiuchus (d^l25 pc) 
and one of the most symmetric, spherical envelopes in 8/im extinction and N2H+ emission, plotted 



in Figure 11 The velocity gradient in this envelope is quite small, illustrated by the channel 
maps with the eastern part of the envelope coming into view just one channel before the rest. 
Figure [12] shows the PV diagram, the emission has a roughly constant linewidth and velocity across 
the envelope and only a slight velocity gradient is evident. Furthermore, the N2H^ emission is 
depressed toward the protostar, as compared to the surrounding emission. This effect appears to 
be a depletion/destruction process since the lines are not optically thick (Paper II). 

This is the only envelope in our sample that can be reasonably described by an axisymmetric 
envelope model; however, it was necessary to lower the mass to 0.1 M©, such that the linewidth 
from infall was not too large. The lower mass is appropriate for this source given its low luminosity 
(L5c»/=0.25 Lq] the envelope does have at least 0.8 Mq of material within 10,000 AU, with increasing 
mass on yet larger scales, suggesting that the central object could be more massive. On the scales 
probed by the CARMA data (R^3125 AU; 2^") the lower-limit on the mass measured from the 
8/im extinction data is just 0.2 Mq. Between 3000 AU and 6000 AU scales, probed mainly by the 
single-dish data from Paper II, the linewidth from spherical infall with dominant envelope mass 
would be ^0.2 km s~^. Thus, the data for this source appear to be consistent with spherical infall 
given the low-mass of the envelope; however, solid-body rotation and filamentary collapse cannot 
be ruled-out by these data. The lack of a significant velocity structure makes the kinematic models 
degenerate for this object. 



3.5. HH211 



HH211 was classified as an irregular envelope in Paper I; however, its structure in N2H^ 



appears filamentary, but not as regular as LI 157 or CB230 (Figure 13). The peak N2II+ emission 
is just southwest of the protostar, coincident with the NH3 peak ( [Tanner fc Arce|[2011 ). The mass 
of this envelope is found to be ^1.1 Mq from 8/im extinction (Paper I) and ~0.65 Mq in the inner 
envelope from N2H+ emission (Paper II). We also find a velocity gradient that is approximately 
normal to the outflow and changing. The PV plot in Figure [T4[ shows that to the northeast, the 
velocity gradient is smaller, while southwest it increases, where is where it blends with an apparent 



second velocity component Paper II and Tanner h Arce (2011). 



On the whole, the data appear most consistent with solid-body rotation, but this is uncertain 
given apparent velocity gradient change and the irregularity of the envelope. The smaller gradient 
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to the northeast is consistent with the infahing filament models, which we overlay in Figure [T4j The 
model is projected by $o = 15° and has Rq — 10 AU.The apparent change of the velocity gradient 
could be due to complex projection effects or the surrounding environment. Even though there is 
substantial mass measured in this envelope with an asymmetric distribution, it is not plausible that 
it could generate such a large linear velocity gradient from infall alone in the context of these simple 
models. Furthermore, the linewidth of the envelope northeast of the protostar is substantially more 
narrow than the linewidth southwest of the protostar. At the N2H^ peak, the linewidth broadens 
substantially and could be related to the outflow as suggested in [Tanner fc Arce (2011), rather 
than inner envelope kinematics. 



4. Discussion 



The direction of the jet and outflow is thought to correspond to the rotation axis of the 
protostar and disk (e.g. Shu et al. 1987). Since circumstellar disks form in the inner envelope 



due to conservation of angular momentum, it has been natural to interpret large-scale velocity 
gradients which are normal to the outflow as rotation (e.g. Goodman et al.|1993 ; Caselli et al.|2002 



Belloche et 'aL||2002t [Belloche fc Andre||2004] [Chen et al.||2007[ and Paper II). If velocity gradients 
probe rotation, they give a measure of the angular momentum in a collapsing envelope, enabling 
the centrifugal radius to be calculated, the radius where a rotationally supported disk will form. 
However, it was shown in Paper I that envelopes are often asymmetric and filamentary and in 
Paper II that the large-scale kinematics are often quite complex as compared to a linear velocity 
gradient. These factors, along with the examples shown in Figure [T] make it uncertain whether or 
not velocity gradients in these envelopes reflect rotation alone. The qualitative arguments shown 
in Figure [l] are reinforced by the simple models of axisymmetric and filamentary collapse shown in 
Figure [2] &|3l 

One of the keys to disentangling the kinematics in protostellar envelopes is higher resolution 
data provided by interferometers. This is because rotation velocities are expected to become most 
important on small-scales due to conservation of angular momentum. Thus, even if the velocity 
structure on larger-scales probed by single-dish observations is not wholly due to rotation, we would 
naively expect rotation velocities to become more important on smaller scales. The interferometric 
PV diagrams and channel maps of these systems reveal small-scale kinematic detail invisible in 
single-dish data, but are also sensitive to emission on size scales between 1000 and 10000 AU. On 
this scale, differences between the various model PV diagrams shown in Figures [2] and [3] become 
apparent. The sources presented in this paper represent some of the most well-ordered systems 
kinematically and morphologically from Paper II. However, in spite of this, there are still difficul- 
ties in disentangling rotation versus infall in the velocity fields and PV diagrams which must be 
recognized when using data of this kind to test theories of star formation. 
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4.1. Velocity Gradients as Rotation 



It was shown in Paper II that the envelopes often have substantial velocity gradients (^2.3 
km s~^ pc~^ ) out to 10000 AU and beyond, measured by single-dish mapping. The majority of the 
envelopes in the sample also had velocity gradients that were within 45° of normal to the outflow 
direction, reinforcing an interpretation as rotation. Previous work has often made the assumption 



that the envelope velocity profile reflects solid-body rotation, even at high-resolution (e.g. Chen 
et al.|200"7 ). The protostars L1165, CB230, and L1157 are inconsistent with solid-body rotation with 
their constant velocities out to large-scales. Therefore, if the line-center velocities reflect rotation, 
then it must be a differential rotation curve with non-constant angular momentum; constant angular 
momentum would yield v oc R~^. HH211 specifically could be consistent with solid-body rotation 
and IRAS 16253-2429 could also be represented by sohd-body rotation but not uniquely. If the 
characteristic large-scale velocity gradient for HH211 from Paper II (6.9 km s~^ pc~^ ) is interpreted 
as rotation, the inferred centrifugal radii of this material is ^2800 AU by applying Equation 5. This 
assumes that material from 10000 AU fell-in to to i?c with a 0.5 Mq central object. However, this 
calculation assumes a constant angular velocity for the entire envelope. 

With the interferometer and single-dish data from Paper II, we can measure the line-center 
velocity from the inner to outer envelope at many positions and calculate the centrifugal radius of 
material throughout the envelope. We have calculated Rc versus radius from the observed line- 
center velocities, assuming a central object mass of O.5M0 and that the observed velocities reflect 
rotation. The results in Figure [15] show that inside of ^5000 AU the centrifugal radii are below 



1000 AU. However, on scales greater than 5000 AU, the centrifugal radii are in excess of 1000 AU, 
except for IRAS 16253-2429. HH211 notably has an unrealistic Rc > 10000 AU at i? - 6500 AU. 
All these centrifugal radii are substantially larger than the characteristic sizes of circumstellar disks 
(-250 AU (l Andrews fc Winiams|[2007| [Andrews et al.||2009D ). We note that Rc could be smaller, 
if the central object masses are larger than assumed. Moreover, if much of the material observed 
on 10000 AU scales does become incorporated into the protostar, then the implied centrifugal radii 
will be smaller by about a factor of 3 at i? = 10000 AU since the masses at i? < 10000 AU — 1 Mq; 
however, this does not fully alleviate the problem since corrections will be less at smaller radii. 

Rotation of the measured magnitude is expected to cause fragmentation during collapse on 
large, -1000 AU, scales (Rafiko^[2005l [2007} [Kratter et al.|[20Tol Zhu et al. 2011). Of the sources 



shown, only CB230 is a wide binary with a separation of —3000 AU. However, the gas kinematics 
around CB230 do not indicate rotationally supported motion on this scale, suggesting that the 
companion did not form via rotational fragmentation. The large centrifugal radii inferred for four 
of five sources, especially HH211, may indicate that another dynamical process is contributing to 
the observed line-center velocities. If the observed velocities are indeed rotation, then to prevent 
fragmentation, the outer envelopes must be removed before they fall-in to the centrifugal radius. 



Outflows are one way to entrain and evacuate ambient material (Arce & Sargent 2006); however 



it seems unlikely that this mechanism can efficiently work for highly filamentary envelopes that 
are extended in the direction normal to the outflow. The physical implications of the observed 
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velocities reflecting only rotation seem to suggest that there are additional contributions to the 
velocity field from dynamical processes other than rotation. 



4.2. Velocity Gradients from Projected Infall 

Since all these systems appear to be gravitationally bound (Paper II), the most likely process 
that could contribute to the line-center velocity field away from the region of outflow influence is 
infall. Figure [l] shows that if an envelope is not axisymmetric, then infall could manifest itself in 
the observed velocity field across the envelope and a velocity gradient normal to the outflow does 
not necessarily imply rotation, also see Figure [3j The models demonstrate that the radial motion 
of a globally infalling filamentary envelope can produce detectable velocity shifts on either side of 
the envelope even with only $o — il5° projection angle within the plane of the sky. The velocity 
shift remains relatively constant on the same scales that we see in our observations. Moreover, if 
the envelope is only marginally resolved (as in the single-dish observations), such a velocity shift 
would likely present itself as a linear gradient in lower-resolution single-dish data. 

CB230, L1157, and L1165 all have velocities that are approximately constant between R^ISOO 
and 10000 AU in the PV plots (Figures [5) [7| & [lO| Table 2). The PV structures agree with the 
predicted kinematic structure for an infalling filament and are inconsistent with what would be 
expected from axisymmetric collapse or solid-body rotation. The resemblance of the PV plots 
to the predictions of the infalling filament model, as well as the filamentary morphology of the 
envelopes, leads us to suggest that the principal dynamic process producing the velocity gradients 
is large-scale infall and not rotation. However, IRAS 16253-2429 could be described by any of 
the kinematic models and HH211 appears most similar to rapid solid-body rotation but the non- 
constant gradient and overly large centrifugal radii may indicate that something more complex is 
taking place. 

The results from the analysis of these data lead us to conclude that many of the large-scale 
velocity gradients observed in Paper II and other studies most likely arise primarily from infall 
velocities being projected along our line of sight, with rotation being a subordinate contribution to 
the velocity field. Furthermore, we suggest that even systems which appear to have linear velocity 
gradients on large-scales probably do not reflect rotation, given the results for HH211. The linear 
velocity gradients probably arise from complex projection effects and the formation of cores in a 
turbulent medium. 



Evidence for large-scale infall has been previously seen in the star-less core L1544 (Tafalla et al 



1998D and possibly NGC 1333 IRAS 4B ( |Di Francesco et al.|2OOl||j0rgensen et al.|2007D . Moreover, 
radiative transfer modeling of the blue-asymmetric line profiles of optically thick tracers indicates 
that infall radii are often of order 5000 AU ( [Walker et aT1|1986| [Zhou et al.||l993| [Ward-Thompson 
& Buckley 2001; Narayanan et al. 2002). Furthermore, starless cores forming in the filamentary 



dark cloud L1517 appear to have material flowing toward the dense cores Hacar &; Tafalla (2011) 
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as viewed in the line-center velocities of the filaments. The large-scale flow of material towards the 
star-less cores embedded within filaments is consistent with our picture of the envelopes infalling 
from large-scales in these protostellar systems. We may be seeing the later evolution of this process 



in our sample, while Hacar &; Tafalla (2011) is seeing the beginning. 



The possibility of projected infall makes the true rotation velocities difficult to disentangle 
in protostellar envelopes with complex morphological structure, limiting the ability to definitively 
measure angular momentum at large-scales. Rotation must be measured at smaller-scales, where 
the rotation velocity will be larger due to conservation of angular momentum during collapse. 



4.3. Effects of Extended Mass 

The CMU infall model assumes a dominant central mass in the calculation of the velocity 
field. However, there is substantial mass on scales larger than ^1000 AU in all sources that could 
certainly affect the kinematics observed. We can simplistically characterize the effects of extended 
mass by employing Gauss's law for gravity, in which any mass interior to it will appear as a point 
mass and the free-fall velocity is 

Since M oc p(i?)R^, we can rewrite the free-fall velocity with the following proportionality 

Vff oc p{Ry/^R. (8) 

Thus, an object with a spherically averaged density profile of p{R) oc will have Vff^ constant 
out to large-scales. Most of the objects examined in this paper have approximately constant 
velocities on large-scales, strongly suggesting that they are infalling and extended mass is influencing 
their kinematics. Even with significant extended mass, the predictions of the filament model remain 
valid; extended mass simply increases the projected infall velocities on large-scales, such that they 
do not approach zero as quickly. In the axisymmetric case, the extended mass would make the infall 
velocities remain quite broad on large scales, inconsistent with the observations. Lastly, projected 
infall could mimic the signature of solid-body rotation if the density were constant with increasing 
radius; however, on the scales of the protostellar envelopes (^0.1 pc) there does not appear to be a 
region of constant density in the molecular line emission or 8/xm extinction. Nevertheless, this effect 
could be evident in larger-scale molecular clouds such as those observed by [Arquilla fc Goldsmith 
(fT986|). 



4.4. Linewidth Constraints on Infall and Morphology 



Evidence for large-scale infall is also found in the linewidths of the envelopes. The N2H+ 
and NH3 linewidths in protostellar envelopes are often observed to be about a factor of two or 
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three larger than expected for purely thermal linewidths (0.13 km s ^ for N2H^ and 0.22 km s ^ 
for NH3) for T = 10 K gas (e.g. ICaselh et al.||2002t |Chen et al.]|2007[ Paper II). Infall motions 
themselves on large-scales, in an axisymmetric envelope, would cause substantial line broadening 
(^0.5 km s~^) due to the superposition of radially infalling material along the line of sight (Figure 
[2]). The line broadening generated by the infalling motions of filamentary envelopes in Figure [s] is 
often just 0.2 km s~^ between radii of 2000 AU and 10000 AU, which is consistent with what is 
observed in the PV plots and in the linewidth fields presented in Paper II. Thus, large-scale infall of 
axisymmetric or non-axisymmetric protostellar envelopes could naturally give rise to the observed 
line broadening and it would not be necessary to invoke a turbulent velocity component to explain 
the observed linewidths. 



4.5. Comparison with Simulations 



Simulations of protostellar core/envelope formation within large-scale molecular clouds also 



support the idea that rotation is not being accurately probed by line-center velocity maps. (Burkert 



fc Bodenheimer|[2000 ) showed that for individual cores, the line-center velocity gradient cannot be 
assumed to probe the angular momentum for individual cores, only the distribution of angular 



momenta can be derived. More recently, Dib et al. (2010) found that angular momenta from cores 



estimated from the line-center velocity field versus 3D velocities are overestimated by a factor of 
^10. However, they only remark that the difference arises from the line-center method fitting 
a global linear gradient and the 3D angular momentum including contributions from "complex 



dynamical behavior," which could be taken to mean turbulence and/or infall. Finally, Smith et al 



(2011) investigated the collapse of filamentary cores in their global simulation, finding that there is 



infall from larger-scales toward the cores themselves. Furthermore, recent detailed analysis of the 
velocity fields show significant complexity with contributions from infall, turbulence, and some net 
rotation (Smith et al. 2011 in preparation). Realistic simulations of collapse have the promise to 
enable a better understanding of the net angular momentum in cores as they will enable statistical 
comparisons to observational data, but with the ability to differentiate rotation and infall. 



5. Summary and Conclusions 



We have presented an analysis of the dynamical processes giving rise to the observed kinematic 
structure in protostellar envelopes. This was done by comparing interferometric N2H+ and NH3 
channel maps and position-velocity diagrams to our qualitative expectations for envelope kine- 
matics and the predicted kinematics of collapse models. We have constructed an analytic model 



that approximates filamentary infall, based on the rotating collapse model (Ulrich 1976; Cassen 



fc Moosman||1981 ), enabling us to explore projection effects in a simple non-axisymmetric system. 
We have matched these models to the observations in order to determine the likelihood of rotation 
or infall in an axisymmetric or filamentary envelope producing the observed velocity structure. 
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The PV diagrams of the infalling axisymmetric and filamentary models are generally distinct from 
each other on scales between 1000 and 10000 AU. The models of filamentary collapse viewed at 
different projection angles show that infall motions along the long- axis of a protostellar envelope 
can produce systematic velocity shifts across the protostar without the need for significant rotation 
at large-scales. These velocity shifts when viewed at low resolution in single-dish measurements 
would give the appearance of an approximately linear velocity gradient. 

The velocity structures of L1157, CB230 and L1165 are found to be inconsistent with uniform 
rotation and tend to have constant velocity outside the inner envelope. If the envelope velocities 
reflected rotation, then the centrifugal radii calculated on scales larger than 5000 AU would be > 
1000 AU, larger than a typical circumstellar disk. The envelope of HH211 has a velocity structure 
consistent with solid-body rotation; however, the envelope is a complex filament and the implied 
centrifugal radii are unrealistically large. IRAS 16253-2429 has the most symmetric, spherical 
envelope in the sample and also has very little velocity structure, making all the kinematic models 
degenerate; however, if rotation is assumed the implied centrifugal radii are reasonable for this 
object. 

Given that the kinematics of LI 157, CB230 and LI 165 seem to reflect projected infall and 
HH211 likely has some infall component to its velocities, we therefore suggest the following scenario. 
The protostellar envelopes we have observed are in a state of global collapse, and the projection of 
the infall velocities along the filamentary envelopes gives rise to the velocity distributions observed 
in the PV plots and moment maps in Paper II. This suggests that envelopes with greater degrees 
of morphological complexity should have similarly complex velocity fields resulting from projected 
infall velocities. 

While there must certainly be net rotation in these envelopes, the degree to which it is solid- 
body is uncertain and global simulations indicate that the turbulent medium from which the cores 
formed would impart net angular momentum to the envelope. While we cannot unequivocally prove 
that our observations are probing infall rather than rotation, we have shown that it is a viable 
interpretation in asymmetric systems. Moreover, the possibility of infall velocity components in 
line-center velocity fields indicates that measurements of angular momentum on large-scales can 
be highly inaccurate. To more definitively measure angular momentum, one must look on small- 
scales where conservation of angular momentum may make the rotation velocities larger and more 
apparent. 
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Sphere/Axisymmetric Filament 




Fig. 1. — These plots schematically show the velocity structure expected to be observed from the 
dynamical processes of rotation, infall, and their combination for particular envelope morphologies. 
The top row demonstrates spherical/axisymmetric envelopes and the bottom row demonstrates 
filamentary envelopes, with schematic position-velocity diagrams below each drawing. The plots 
illustrate that solid-body rotation in an axisymmetric envelope {top left) will show a linear velocity 
gradient and infall alone {top center) only increases linewidth. The convolution of these processes 
is shown in the top right panel. A filament also shows a linear velocity gradient for uniform rotation 
{bottom left); however, infall along a filament turned toward our line of sight would generate an 
approximately linear velocity structure in the absence of rotation {bottom middle). When these 
process are convolved in the bottom right panel, the two processes are difficult to separate. 
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Axisymmetric Axisymmetric Axisymmetric Rc = 0100 
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Fig. 2. — Plots of infall streamlines and projected velocities for an axisymmetric envelope with 
Rc=WO AU. The left panel shows the streamlines out to 10000 AU while the middle panel zooms 
in on the inner 1000 AU. The circle at the center represents the edge of the circumstellar disk 
forming at Re- The right panel shows the predicted PV structure for the axisymmetric envelope. 
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Fig. 3. — Plots of infall streamlines for a filamentary envelope with Rc=lOO AU. From top to 
bottom, filamentary envelopes are shown with projection angles of $0= 0, -30, and 30°. The left 
panel shows the streamlines out to 10000 AU while the middle panel zooms in on the inner 1000 
AU. The circle at the center represents the edge of the circumstellar disk forming at Rc- The 
right panel show the PV structure of the filamentary streamlines. The filamentary envelopes at a 
particular rotation angle can be thought of as showing a selected portion of the axisymmetric PV 
diagram. 
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Fig. 4.— L1157- Channel maps of the isolated N2H+ {JFiF = 101 012) line; the blue and 
red arrows mark the outflow direction and the cross marks the location of the protostar. The 
higher velocity data in the top and bottom rows mainly trace the likely outflow entrained material, 
while the middle two rows primarily reflect emission due to envelope kinematics. Note the velocity 
gradient of the large-scale material from 2.45 km s~^ to 3.075 km s~^. The large-scale emission 
on the east and west sides of the envelope has a rather constant velocity. The extension of N2H+ 
emission to the south traces an outflow cavity wall whose kinematics also reflect some outflow 
entrainment. Contours in the channel maps start at ±3cr increase by ±3cr intervals, where a=0.08 
K. The contours in the integrated intensity map are ±6, 9, 15, 30, 60a, where a=0.02 K km s~^. 
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Fig. 5. — L1157- The left panel shows the IRAC 8/im image with PdBI N2H+ integrated intensity 
contours overlaid of the isolated N2H+ line {JFiF = 101 012). The dashed lines mark the regions 
where the position-velocity cut was taken and point to respective ends of the PV plot in the right 
panel. The position of the protostar/continuum source is marked with a white cross. The position- 
velocity cut is shown for the isolated N2H+ line, showing that the two lumps of N2H+ emission are 
directly associated with regions of increased linewidth, likely due to outflow effects. Outside of this 
region, the linewidth becomes fairly narrow in the outer envelope with constant velocity. The green 
lines overlaid are an filamentary collapse model with i^^^lOO AU and a projection of $0 = 15°. 
The contours in the integrated intensity map are ±6, 9, 15, 30, 60cr, where cr=0.02 K km s~^; 
the contours in the PV diagram are zbScr and in increments of 3a where (j=0.17 K. 
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Fig. 6.— LI 165- Channel maps of the isolated N2H+ {JFiF = 122 Oil) line; the blue and red 
arrows mark the outflow direction and the cross marks the location of the protostar. There is a 
rather well-behaved global velocity structure in LI 165 from blue-shifted northwest of the protostar 
and red-shifted to the southeast. There does not appear to be significant outflow effects on the 
kinematic structure. There is also higher velocity red-shifted emission apparent on small-scales 
near the protostar, but very little blue-shifted emission with higher velocities. The contours in the 
integrated intensity map are ±2.5, 3, 6a, where cr=0.207 K km s~^; the contours in the channel 
maps are ±3, 6a, where a=0.4 K. 
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Fig. 7.— LI 165- The left panel shows the IRAC 8 jam image with CARMA N2H+ integrated 
intensity contours overlaid. The dashed hnes mark the regions where the posit ion- velocity cut 
was taken and point to respective ends of the PV plot in the right panel showing the N2H+ 
(JFiF = 122 Oil) transition. The position of the protostar/continuum source is marked with a 
white cross. Northwest of the protostar the lines tend to be more narrow than on the southeast side 
of the protostar where the N2H+ emission is both broad and red-shifted. The green lines overlaid 
are an filamentary collapse model with Rq^W AU and a projection of $0 — 15°. The contours in 
the integrated intensity map are ±2.5, 3, 6cr, where cr=0.207 K km s~^; the contours in the PV 
diagram are ±3cr, increasing in increments of 3a where a=0.225 K. 
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Fig. 8.— L1165- The left panel shows the IRAC 8/xm image with CARMA HCO+ (J=l^ 0) blue 
and red-shifted emission, summed over -3.5 to -2 km s~^ and -1.0 to km s~^, plotted as blue and 
red contours respectively. The contours levels are 3, 6, and 8.25cr ((j=0.175 K) for the blue-shifted 
emission and 3, 9, and 18a (a=0.212 K) for the red-shifted emission. The blue and red-shifted 
emission from HCO^ is located symmetrically about the protostar, normal to the outflow. The 
dashed lines mark the regions where the position-velocity cut was taken and point to respective 
ends of the PV plot in the right panel. The position of the protostar/continuum source is marked 
with a white cross. The position-velocity cut shows that the blue and red-shifted emission traces 
higher velocity material and there is a slight gradient of material going to higher velocity closer to 
the continuum source. The line center velocity of -1.5 km s~^ is plotted as the dotted black line. 
The green lines overlaid are an filamentary collapse model with Rc^W AU and a projection of $o 
= 15°. The PV plot contours start at 3a and increase in 3a intervals (cr=0.2). 
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Fig. 9. — CB230- NH3 integrated intensity map of the main NH3 (1,1) lines {upper left) and 
channel maps of the two main NH3 (1,1) lines are shown in the rest of the panels. The protostars 
are marked with crosses and the outflow direction is drawn from the main protostar. The emission 
between 2.54 and 2.85 km s~^ consists of blended emission from both lines, while the remaining 
emission is mostly free of blending. The channel maps show a clear velocity gradient across the 
source, appearing normal to the outflow. Notice how the emission tends to avoid the location of 
the main protostar, owing to destruction of NH3 near the protostar. The contours in the integrated 
intensity map start at ±3a, increasing in increments of 3a and a=0.195 K km s~^; the contours in 
the integrated intensity map are ±3, 6cr, where a=0.325 K. 
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Fig. 10. — CB230- Same as Figure [s] but with VLA NH3 (1,1) observations and the secondary 
source is also marked with a white cross. The PV diagram in the right panel appears different from 
those constructed from N2H+ because the NH3 emission consists of two blended hyperfine lines; the 
NH3 satellite lines with greater separation are shown in the PV plot rather than the more blended 
main lines. The line-center velocity from the NH3 emission traces an abrupt velocity shift from 
red to blue-shifted emission. The transition region starts at the protostar and does not finish until 
+10'^ The NH3 linewidth does not show much detail other than having is peak coincident with the 
highest intensity NH3 emission. The green lines overlaid are an filamentary collapse model with 
Rc^W AU and a projection of $0 — 20°. The contours in the integrated intensity map start at 
zbScr, increasing in increments of 3a and a— 0.195 K km s~^; the contours in the PV diagram start 
at ±3cr, increasing in increments of 3a and cr=0.11 K. 
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Fig. 11.— IRAS 16253-2429- Channel maps of the isolated N2H+ (JFiF = 122 Oil) line; the 
blue and red arrows mark the outflow direction and the cross marks the location of the protostar. 
IRAS 16253-2429 is the envelope that is most nearly round. The velocity gradient across the 
envelope is very small, but is clearing seen in the channel maps with emission on the eastern side 
of the envelope coming into view at blueshifted velocities and the western side with redshifted 
velocities. The green lines overlaid are an axisymmetric collapse model with Rc^lOO AU and Mq 
= 0.1 Mq. Contours in the channel maps start at ±3cr increase by ±3cr intervals, where a=0.08 
K. The contours in the integrated intensity map are ±2.5, 3, 6cr, where cr=0.207 K km s~^. 
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Fig. 12.— IRAS 16253-2429- Same as Figure [s] but with CARMA N2H+ observations of the 
{JFiF = 122 Oil) transition. The N2H+ emission matches up quite well with the 8/xm extinc- 
tion. Furthermore, the N2H+ appears peaked around the protostar and not coincident with it. This 
is shown as the decreased emission in the PV plot in the right panel. The velocity gradient in the 
PV plot is minute but is present. The green lines overlaid are an filamentary collapse model with 
Rc^W AU and a projection of $0 — 15°. The contours in the integrated intensity map start at 
zb3(j, increasing in increments of 3a and a— 0.22 K km s~^; the contours in the PV diagram start 
at zb3cr, increasing in increments of 3cr and cr=0.25 K. 
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Fig. 13.— HH211- Channel maps of the isolated N2H+ {JFiF = 101 012) line; the blue 
and red arrows mark the outflow direction and the cross marks the location of the protostar. The 
N2H+ emission in the different velocity channels clearly trace a filamentary envelope with a velocity 
gradient across the protostar that appears fairly linear. Contours in the channel maps start at zh3cr 
increase by =b3cr intervals, where cr=0.4 K. The contours in the integrated intensity map are ±10, 
20cr, where cr=0.149 K km s"^. 
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Fig. 14.— HH211- Same as Figure [s] but with CARMA N2H+ observations. The PV plot in the 
right panel is constructed from the N2H+ {JFiF = 101 012) transition. Notice how the slope 
of emission in the PV diagram changes from one side of the protostar to the other. The green 
lines overlaid are an filamentary collapse model with Rc^lO AU and a projection of $0 = 15° that 
are meant to represent the kinematics on the northeast side of the envelope. The contours in the 
integrated intensity map start at zblOa, increasing in increments of lOcr and cr=0.149 K km s~^; 
the contours in the PV diagram start at ±3a, increasing in increments of 3a and a=0.2 K. 
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Fig. 15.— Plot of centrifugal radius (Re) versus radius for L1157, L1165, CB230, IRAS 16253-2429, 
and HH211 assuming the observed velocities reflect rotation. Rc is calculated from the observed 
line-center velocities from Paper II and converted to a centrifugal radius using the relation Rc = 
R'^{V - Visr)'^/GM (Equation 4), where M = 0.5 Mq is assumed. 
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Table 2. Model Parameters 



Source 


Model 


Mcent 


Rc 


Proj 


. Angle 






Mq 


(AU) 




(°) 


L1157 


Filament 


0.5 


100 


15 




L1165 


Filament 


0.5 


10 


15 




CB230 


Filament 


0.5 


10 


20 




IRAS 16253-2429 


Filament / Axisymmetric 


0.1 


100 







HH211 


Filament 


0.5 


10 


15 





Note. — 



Parameters of models which describe the PV data. 



